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We study a hierarchy of directed percolation (DP) processes for particle species A,B,.. ., unidirec- 
tionally coupled via the reactions A — > _B, . . .. When the DP critical points at all levels coincide, 
multicritical behavior emerges, with density exponents which are markedly reduced at each 
hierarchy level k > 2. We compute the fluctuation corrections to to 0(e = 4 — d) using field- 
theoretic renormalization group techniques. Monte Carlo simulations are employed to determine 
the new exponents in dimensions d < 3. 
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Nonequilibrium critical phenomena have been the fo- 
cus of intense theoretical studies in recent years, promi- 
nent examples being phase transitions in driven diffu- 
sive systems [|l|, kinetic Ising models nonequilibrium 
growth models Q, diffusion- limited chemical reactions 
percolation-like processes ||], and generally in lat- 
tice models or cellular automata. The concept of only 
a few distinct universality classes which determine the 
critical exponents has been remarkably successful for the 
description of static and even dynamic critical phenom- 
ena near equilibrium phase transitions |^]. However, it is 
still an unresolved issue whether a similarly small number 
of universality classes can be identified in nonequilibrium 
situations, which potentially contain much richer behav- 
ior. In addition to the dynamic universality classes listed 
in Ref. 0, some notable new candidates have emerged 
in quite different circumstances, for example, the rough- 



levels, respectively. 



, and the critical 
as described by 



ening transition in the KPZ equation 
point of directed percolation (DP) | 
Reggeon field theory jl^]. In fact, the DP universality 
class appears to cover the majority of phase transitions 
from non-trivial active into absorbing states where the 
order parameter noise vanishes. However one striking 
exception to this rule occurs when the relevant dynamic 
processes are constrained by an additional local "parity" 
symmetry, in which case the universality class is that of 
branching and annihilating random walks with an even 
number of offspring ||ll|,|l^ . 

This work is originally motivated by recent studies of 
a nonequilibrium growth model for adsorption and des- 
orption of particles which displays a roughening transi- 
tion in d = 1 ||l^. The key feature is that desorption 
may take place only at the edge of an existing plateau 
so that particles cannot be desorbed from a completed 
layer. Because of this property the dynamic processes 
at a given height level decouple from all higher levels, 
i.e. the processes at different heights are unidirection- 
ally coupled. This growth model can in fact be related 
to a reaction-diffusion problem for a hierarchy of parti- 
cle species A,B,..., which correspond to different height 



Each of these species is subject to 



the prototypical DP processes (see Refs. |10 
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with reaction rates cta, Xa, and ^a, respectively (and 
accordingly for particles B, . . .). However, there is also a 
coupling from each previous hierarchy level to the next 
one via random particle transmutations 



A^ B , B 



C 
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with rates /.iab, fJ-BCi ■ ■ but no feedback mechanism 
from the B to the A species etc. For simplicity, we 
shall assume that the diffusion constant D for all par- 
ticle species A, B, . . . is identical. When ^ab = Mbc = 
. . . = 0, the system decouples, and for the i-th species 
there is a critical point separating an active from an ab- 
sorbing state, where essentially the branching rate at and 
decay rate fii balance one another. At this point there 
are three independent critical exponents governing (i) the 



divergence of the correlation length ^ oc 



(where 



Ti cx fii — ai denotes the deviation from the critical point), 
(ii) the critical slowing down of characteristic time scales 
tc cx oc jril""!!, with — ziy±, and (m) the value 
of the asymptotic particle density in the active phase, 
ni{t — > CX)) cx l^il^. Alternatively, this third independent 
exponent may be swapped for r]±, which characterizes 
how the equal-time pair correlation function decays at 
the critical point n = 0, G(|x|) cx l/|x|^+^-2-i-»u [gee 
Eq. ( pO| ) below] . These exponents should be those of the 
DP universality class, with upper critical (spatial) di- 
mension dc = 4. Numerical simulations in d = 1 have 
confirmed this, and furthermore revealed that i^^ ~ 1.1 
and z « 1.6 remain unchanged at all hierarchy levels 
even when the rates fiAB ■ ■ ■ are switched on. However, 
while Pa = P^^^ ~ 0.27 as in DP for the primary A 
particles, the particle density exponents for secondary 
and higher hierarchy levels are considerably reduced to 
[3b = » 0.11, and /J^^) « 0.04, provided the critical 
points for the processes on the different levels coincide 
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d= 1 


d = 2 


d = 3 


d = 4- e 




0.271(10) 
0.108(10) 
0.038(8) 


0.57(5) 
0.30(4) 
0.13(3) 


0.78(7) 
0.35(6) 
0.15(4) 


1 - e/6 + 0(e') 
l/2~e/Q + 0{e^) 
l/4-e/6 + 0(e2) 



TABLE I. Results for the particle density exponents /S'-'''' 
obtained from Monte Carlo simulations and RG calculations. 



In the following, we shall mainly study the two-level hi- 
erarchy of DP processes (|l|) for A and B particles coupled 
by the reaction A ^ B. We explain the reduced value of 
Z?'-^-' as a consequence of the multicritical behavior that 
emerges when r^i = = r ^ 0. It appears then rather 
likely that similar multicritical behavior should emerge 
in more general contexts than the specific growth model 
and diffusion- limited chemical reactions (0), (H); namely 
whenever DP-likc processes are coupled unidirectionally 
without feedback. A qualitatively correct description is 
obtained by analyzing the corresponding mean-field rate 
equations, which yield [3'^^^ = l/2^~^ at the multicritical 
point on the fc-th hierarchy level. With the aid of the 
renormalization group (RG), we then identify an addi- 
tional independent crossover exponent (p related to the 
relevant scaling field {^ab /D) of the two-level hierarchy 
multicritical point. A scaling relation will be derived ex- 
pressing in terms of (j) and the conventional DP ex- 
ponent (3 — (3^^\ Using a field-theoretic representation 
(see Refs. p^ , p^ ) based on the master equation for this 
stochastic process, we furthermore compute the fluctua- 
tion corrections to and to first order in e = 4 — c? 
1^^ . Finally, these exponents are determined accurately 
in c? < 3 dimensions by means of Monte Carlo simula- 
tions. 

We start by writing down the mean-field rate equations 
for the average local densities ri/i(x, t) and ^^(x, i) of the 
A and B particles, both subject to the reactions (|^), and 
coupled via the random particle transmutation A ^ B, 

dtnA= D{V^ - rA)nA - A^^i , (3) 
dtHB^ D{V^ - rB)nB - \Bn% + i^AsnA , (4) 

where ta = (ma + I^ab - crA)lD and tb = - o-b)ID. 
Considering first the equation for ua, it is clear that for 
> the only stationary state is ua = (inactive 
phase). On the other hand, if r^i < the density will 
asymptotically approach n'^ = D{—rA)/^A > (active 
state), and hence /3^^^ = 1 in this approximation. From 
the diffusive character of Eq. (||) , and also from the func- 
tional form of its linear term, it follows that z = 2 and 
u± = 1/2 in mean- field. These last two results remain 
intact for the B particles as well, even in the presence 
of particle influx via the reaction A B. However, the 
possible stationary states for ub actually depend on the 
value of rA as well as on vb, if ^J-AB > 0. If > 0, 
and therefore n^(i — > oo) = 0, Eq. (^) reduces to a 
mean- field DP process with critical point vb = and 
(3 — 1. On the other hand, in the A-particle active 



state, rA < 0, the stationary solution of becomes 
= [iDrB/2XB)^ + liABu'f IXb^'^ - DrB/2\B\ i.e., 
there is merely a crossover from a low-density to a high- 
density regime in the vicinity of = 0. An interesting 
situation arises when rA = rs = r, as in this case the 
second term in the square brackets dominates for |r| — > 0, 

n'S = [D^iAB{~r)/XAXB]'^^ + 0{\r\) . (5) 

Hence /J^^-' = 1/2 instead of the mean-field DP value. 
More generally, this multicritical behavior arises in the 
regime {DrB/2\BY < D{-rA)fiAB/>^A>^B for T 0. 
On the other hand, for > and [DrB /2\bY ^ 
D{-rA)iJ.AB/><A>^B one finds n"^ = {-rA)tJ'AB/rB\A, 
and thus one expects DP transitions for the B species 
both as rs ^ (with > 0) and ^ (rs > 0) [0. 

Notice that there are three critical lines, namely (r^ = 
0, rB > 0), {rA = 0, rs < 0), and (r^ > 0, = 0), two 
of which are critical lines for the B particles, meeting at 
the point rA — rs — 0. This special point can therefore 
be interpreted as a multicritical point, with the mean- 
field order parameter exponent halved as compared to the 
ordinary critical point ||l^ . Clearly the generalization of 
Eqs. (|), (!) to k hierarchy levels leads to /J^*^) = l/2'=^i 
at the multicritical point = = . . . = 0. 

For d < dc — A, one expects that these mean-field 
values will be strongly modified by fluctuation effects. 
We have performed Monte Carlo simulations in d < 3 
dimensions in order to determine the survival probability 
exponent /J^'^^ /'^|'^^ ^^"^ third hierarchy level. The 

DP processes dJ) are realized by a cellular automaton 
which evolves by parallel updates in which a site at time 
t + 1 becomes active with probability p provided that 
the same site, or at least one of its nearest neighbors, 
was active at time t. The coupling (^) between levels is 
incorporated by the rule that active sites at a given level 
simultaneously impose active sites on all higher levels. 
We perform dynamic Monte Carlo simulations ||l^ where 
one measures the temporal evolution of the system at 
criticality starting from a single seed (a single A particle). 
Assuming that = is identical on all levels k, we 
can use the survival probability to determine the density 
exponents, giving in d = 1 P^^^^ = 0.271(10) as in DP 
!§, = 0.108(10), and /J^^) = 0.038(8). These results 
are consistent with direct estimations of /3^'^) obtained in 
off-critical simulations fl^ . The mean-square spreading 
exponents yield the exponent z, with the results z^^"^ = 
1.57(2) (DP), ^ 1.58(2) and z^^) = 1.59(2), which 
confirms that these exponents are identical on all levels. 
Similar results are obtained in d = 2, d = 3 (see Table I). 

In order to better understand the very low values found 
for P^'''> in the simulations, we clearly have to take fluctu- 
ation effects into account. In low dimensions, where re- 
action noise in the processes (|l|), (||) is highly important, 
a straightforward Langevin-type description of the noise 
can be dangerous, because of the highly non-trivial form 
of the noise correlations. Instead, it is useful to start 
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from the corresponding master equation, and utiHze a 
standard formahsm involving a second-quantized bosonic 
operator representation to derive an effective field theory 
directly It is important to note that apart from 

the continuum limit, this procedure is exact and requires 
no additional assumptions regarding the precise form of 
the noise. In the present case, the field theory describing 
the two-level coupled reactions in d dimensions is defined 
by the action (we omit terms related to the initial state) 

5= j d'^x j dt^d [dt + D{rA - V^)] a - UAd^a 
+\a {da^ + a^a^) — HAsba (6) 
+b [dt + D{rB - V^)] b - aBb% + Xb {bb^ + bH'') } . 

We remark that (a) = ua and (6) = n^, whereas the 
fields a, b have no direct physical interpretation. 

It is now convenient to rescale the fields accord- 
ing to a = (AA/crA)^/^'0o, a = {aA/ XaY^'^iPo, b = 
(A_B/crB)^/^(^o, b = [ob/ XbY^'^^0, and also define new 
couphngs Mo = 2(crAA^)^/^, Uq = 2((j_bAb)^/^, as well 
as /io = ^Asio'A^B / c^b^aY^^ (henceforth, the subscript 
"0" denotes unrenormalized quantities). If we introduce 
a length scale and correspondingly measure times 
in units of (i.e., [Dg] = k°), we find that the new 
fields have scaling dimension k^/^, while [ta] — [tb] = 
[^o] = which are thus relevant perturbations in the 
RG sense. On the other hand, [uq] = [mq] = , and 

the corresponding DP nonlinearities become marginal in 
dc = 4: dimensions, as expected It is important to 
note, however, that [A^] = [As] = k^^*^, and hence these 
couplings are irrelevant as compared to uq and Uq, and 
may be omitted in the effective action (see Ref. [p^), 
which consequently reads 

5eff= fd^'x jdt{i,o [dt + Do{rA ~~ V^)] Vo 

(7) 



u' } 
+ifo [dt + Do{rB ~ V^)] 'Po- Y i'^o'Po - (fofl) | • 

The saddle-point equations corresponding to this action 
are precisely the mean- field rate equations (^), (Q) but 
with ipQ, (fo instead of riA, riB, and with fiAB, ^A, and Xb 
replaced by /i, uo/2, and Mo/2. We remark that the har- 
monic part of this action can be diagonalized only when 
fA 7^ ^B, leading to new cubic vertices mixing the fields 
■00 and (fiQ [compare Eq. (^) below]. However, these ad- 
ditional vertices have no influence on the behavior of the 
A and B species |2^, and the diagonalized theory leads 
to the anticipated DP critical behavior at both critical 
lines (rA = 0, rs 0) and (r^ > 0, = 0) 0. 

Here we are predominantly interested in the multicrit- 
ical point rA — rB = rg ^ 0, and thus we must work 
with the non-diagonal action (|^). Furthermore, due to 
the inclusion of the relevant scaling field fiQ, a number 



of "mixed" cubic vertices are additionally generated at 
the tree level, all of which have scaling dimension k^~'^/^. 
These vertices must be taken into account separately, and 
hence we need to replace (|^) with Smc = S'cff + AS', where 



d'^a 



dt 



(8) 



It is now a straightforward task to compute the renor- 
malized couplings and the RG fixed points to one-loop 
order. The DP nonlinearities in (^ remain unaffected by 
the "mixed" vertices (H), and the stable nontrivial fixed 
points for the associated dimensionless renormalized cou- 
plings U = UoAy^K('^-4)/2 g^^^ ^, ^ u(,v4y\('*-4)/2^ 

where Ad = r(3 - d/2) , are 



[{u/Drf^[{^,'|Drf^Ae|i + 0{^) 



(9) 



Hence the critical exponents 77_l = — e/12 -I- O(e^), i/j^^ = 
2-e/4-f 0(e2), and z = 2 - e/12 -h ©(e^) remain those of 
the DP universality class at all hierarchy levels, whereas 
the exponent /3^*''-' is unchanged only at the first level 

/?(!) = !^{d + z-2 + r/^) = 1 - e/6 + ©(e^) . (10) 

For the similarly defined renormalized "mixed" cubic ver- 
tices, one finds two RG fixed lines. The first is given by 

{s/Or - -{S'/DT, {S/DY [s'/DT = 2{slD)\ (11) 
[{s/DYf = 2{elZf'\slDr + {s'/DT] , (12) 

with stability matrix eigenvalues 0, 0, — e/3, — e/3, which 
imply that this fixed line (including the Gaussian fixed 
point for AS") is unstable for d < 4. The second fixed 
line, 

(.7^)* = {S/DT, {s/DT + [s'/DY = 2(6/3)l/^ (13) 

with again Eq. ( |l^ ) , turns out to be stable as its stability 
matrix has eigenvalues 0, e/3, e/3, 4e/3 |17|. We remark 
that both fixed lines ( pi] ) and ( [l^ ) with (|12|) satisfy the 
condition {e/iY/\s'/D)* + {3/D)*] = -{s/ D)* {s' / D)\ 
which ensures the cancellation of strongly singular (UV- 
divergent in d = 2) diagrams for the renormalization of 
/io. The RG eigenvalue of at these fixed lines then 
becomes — 2 -f- e/6, and j/^ = 2 — e/6, respectively. 
In principle a product of quartic vertices and fiQ might 
enter the renormalizations of the three-point functions 
as well. However, we have checked that these additional 
couplings all have negative RG eigenvalues and are there- 
fore irrelevant [ pTf . 

We next define the crossover exponent cj) related to the 
new scaling field /i/D, as it appears in the general scaling 
form for the B-species "density" field. 



^{\rlf,/D) = \rf' 0i\r\-^^^/D) 



(14) 
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and identify 



(15) 



at the stable fixed line. 

Finally, we relate the above exponent (p to the density 
exponent on the second hierarchy level. We first use 
the RG flow equations for the scale-dependent renormal- 
ized parameters, and then match to our earlier mean-field 
result (|^) for the density in the active phase in the vicin- 
ity of the multicritical point. In terms of the rescaled 
fields, we found = {2hqiI)^ /u'^f/^ for ro ^ 0. Es- 
sentially, we now have to replace the bare parameters 
here by their flowing renormalized counterparts, when 
inverse length scales are rescaled according to k ^ ni. 
To that end, we first note that the fields ip and scale 
as £('i+2~2+')i)/2 [which, with the matching condition 
I = (— r)''^, actually implies the scaling relation cited in 
Eq. (p^]. Furthermore, the 0(e) fixed points require 

\t)/D{l) const., and according to our definitions, 
we have rit) rf-i/"^ and n{i)/D{i) oc e-'*'/"^. Upon 
identifying i = (—r)'^^, we finally arrive at 

p(2) ^ _ = 1/2 - e/6 + 0(e2) , (16) 

using our above results |l^] . Equivalently, this follows by 
demanding that the scaling function in Eq. (^4|) behave 
as 'f'{x) x^/'^ for x — > oo, as prescribed by mean-field 
theory d). The 0(e) result (|l|) constitutes a sizable 
downward renormalization through fluctuation effects, as 
is also evident from the simulation data in Table I. These 
can be explained with the values (f) ~ 0.33, 0.54, and 0.86 
in d = 1, 2, and 3, respectively. 

In principle, these considerations are readily general- 
ized to further hierarchy levels For example, al- 
though a combination of processes A B and B ^ C 
will generate the reaction A — » C as well, one finds that 
the additional relevant scaling field fiAC does not en- 
ter the leading contribution to the C density near the 
multicritical point ta = tb — rc — r —> 0. Hence 
nc ~ [D^^AB^J■%ci-'^)/^A^B>^cV^'^, implying that no 
additional independent exponents need to be introduced. 
We then find /S^^) = _ 30/4 = 1/4 - e/6 + 0{e^). 

In summary, we have studied a hierarchy of unidirec- 
tionally coupled DP processes (no feedback), as defined 
by the diffusion- limited reactions (|l|) , (^) . Already within 
mean-field theory, we have identified a purely dynamic 
multicritical point, which occurs when the DP critical 
points at each level coincide. Using a field-theoretic rep- 
resentation, we have evaluated to 0(e) the crossover ex- 
ponent related to the relevant scaling field ^ab/D. We 
then derived a scaling relation which allowed us to deter- 
mine In accord with our Monte Carlo simulation re- 
sults, fluctuation effects lead to a strong downward renor- 
malization of the density exponents as compared with the 
mean-field values /3''''-' ~ 1/2''^-^. We emphasize that the 
existence of such a multicritical regime is to be expected 
generically whenever processes of the ubiquitous DP uni- 
versality class are coupled unidirectionally. 
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